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We show that a stable spontaneous current appears in a system of a spin density wave 
order with merons, vortices in the spin configuration with winding number +f . A meron in 
the spin density wave order modifies the boundary condition for eigenfunctions around it to 
a sign-change one. As a consequence, two types of stable current states, one with a clockwise 
circulation and the other with a counterclockwise one, arise. When a magnetic field is present, 
one produces a diamagnetic current is chosen. A collection of such currents results in a large 
diamagnetic current; and if the meron density is sufficiently large, a perfect diamagnetism 
is realized. The stability of this diamagnetic current is attributed to the topological nature 
of the merons, and as long as the distribution of the merons remains the same the current 
will persist even in a macroscopic system. 
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It is known that a quantum number for a circular motion is modified from integers to half- 
odd integers when the boundary condition on eigenfunctions is altered from a single-valued 
one to sign-change one by a Berry phase. ^ Such a change was first found in the context of 
the dynamical Jahn- Teller effect,^' ^ but now similar effects have been found in many other 
systems, including band twisting ones.^'^ One of the important consequences brought about 
by this quantum number modification is that the ground state becomes degenerate.^' ''"^ 

To see this, let us consider a particle in a ring-shaped system given in Fig. 1. We denote 
the polar angle for the position in the ring by (f). When the ordinary single-valued boundary 
condition is imposed, eigenfunctions for the particle state are (27r)~^/^e*"'^ {n is an integer). 
The ground state is the n = state, which is not degenerate. On the other hand, eigenfunctions 
for the sign-change boundary condition have the same functional form but with half odd 
integers for n; thus, the ground states are those with n = ib^ and degenerate. If the particle 
has an electric charge, this degeneracy implies an electric current generation in the ground 
state. 

In the present work, we deal with a spin density wave (SDW) ordered state with merons; 
and examine the appearance of circular currents around them due to the sign-change boundary 
condition, a mechanism similar to the one explained above. 

Merons are vortices in the spin configuration with winding number (if the winding 
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Fig. 1. A ring-shaped system and the definition of (f>. 



number is —1 it is cahed, "antimerons").^'^ A typical meron in the xy plane is given by 

^(r)=tan-i|^, (1) 

where {xc, He) is the center of the meron, and rj is the angle for spin direction in the xy plane. 

There is compelling evidence that merons exist in the background antiferromagnetic spin 
order in the pseudogap phase of cuprate superconductors.^'^ If we approximate this background 
spin order as an SDW order, the situation is represented as a system with an SDW order 
with merons. Furthermore, in the pseudogap phase, vortex-like excitations and enhanced 
diamagnetism have been observed; they may be explain by currents around merons we will 
discuss in the present work. 

Now let us consider the ground state vector for an SDW ordered state given by 

|$SDw(0)=n(^°s^ke-*^^'^kT+sin^ke'«/'4H<3|)|vac), (2) 

k 

where cj^^ is a creation operator for a state with wave vector k and spin cr, Q is the Fermi 
surface nesting vector, and is the phase of the SDW order parameter. 

A salient feature of this state vector is that two states with different ^ values are physically 
distinct unless the difference is a multiple of 2tt. Thus, the ground state is degenerate and 
labeled by ^. This degeneracy gives rise to a possibility to have states with a spatially and 
temporally varying ^. 

Using a fermion field operator 

^<T(r) = ^Cka0ka(r), (3) 
k,(T 

where basis functions are given by (pkai^) = V^^^'^e'^^'^ {V is the volume of the system), the 
spin density for |$sdw) is calculated as 

s^{r) = ($sDw|^ (V^jW^/'iW +V'|(r)t/;|(r)) |$sDw) 

= |s|cos(Q • r + ^), (4) 

Sy{r) = ($sDw|^ (V^}(r)v;^(r) - ^/;J(r)^|(r)) |$sDw) 

= |s|sin(QT + 0, (5) 

s,{r) = ($sDw|^ (V'{(r)?/'T(r) - V'lW^/^iW) I^SDw) 
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= 0, (6) 

with |s| given by 

|s| = y-^^cos^ksin^k- (7) 

k 

For spatially modulating ^, the state vector Eq. (2) is modified as 

|5*SDw) = n(^°^^k(r)4-r + sin6'k(r)4:+Qi)|vac), (8) 

k 

where r is the position operator. It is replaced by the argument of a field operator 

with which the contraction of the accompanying creation operator is taken; for example, 

. I 1 ^ I 1 

Vi(r)cos6'k(f)4t = V'(r)cos6'k(r)4-f• 
New fermion operators cj^^ and cj^j^ are creation operators for basis functions that take 
into account the spatial variation of ^, 

^kT(r)=e^5(x{r)-«(r))^^^(r) (9) 

and 

^ki(r) = e*^(>^W+««)0ki(r), (10) 
respectively; and the field operator is now given by 

= XI ^ka(r)cka, (11) 
k,cr 

where the U{1) phase factor e'^^"")/^ that does not alter the spin density is added. Without it, 
the new basis functions are not single- valued: because the value of the phase ^(r) changes by 
27r after a circular transport around a meron, the sign change of e*^^"")/^ and e"*^^'")/^ occurs; 
thus, e*^'^'')/^ is needed to compensate this sign change. The compensation may be done, for 
example, by choosing % to be x(r) = ^(r), but the real % solution should be determined by 
the minimum total energy condition. Note that if xi is a solution, —xi is also a solution by 
the time-reversal invariance of the Hamiltonian. 

The appearance of the phase factor e^^^^^^^ may also be understood in the following way: 
instead of Eq. (8), the straightforward inclusion of the spatially varying ^ phase yields 

I^SDw) = n*^'^°^^k(r)CkT +s™^k(r)4+Qi)|vac), (12) 

k 

with fermion operators cj^.^ and c]^^ for 

^kT(r) = e-i«W)0kT(r) (13) 

and 

0ki(r)=e^l^«<^ki(r), (14) 
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respectively. Then, the wave function 

t(ri,--- ,rjv;^(r)) = (ri,--- ,rjv|^SDw) (15) 

is a multiple- valued one due to the multi-valuedness of e*^^'")/^. 

The total wave function, a product of a wave function for the above internal state and a 
wave function for a motion of a system as a whole, is given by 

e^E,x(r,)^(ri,---,r;v;e(r)), (16) 

where e^^^^^""^^ is the wave function for the whole system motion. The total wave function 
must be single- valued; thus, x must have the multi-valuedness that compensates the muti- 
valuedness of ^. 

If we choose x to be |x, the single- valued total wave function is obtained; and the resulting 
wave function is, 

eiE,|x(r,)^(ri,--- ,rjv;.e(r)) = (ri,--- ,r^|$SDw)- (17) 
Therefore, e*^^ 2><:(rj) ^.^^^ -j^^ regarded as a wave function for the motion of the system as a 
whole. 

The nontrivial phase factor e^'X-^^^l'^ gives rise to a circular current around the meron. We 
will call it, a "meron current". Corresponding to the two solutions x\ ^^id — Xi, there are two 
choices in the direction of the current. 

In the following we assume that the spatial variation of ^k's are negligibly small. 

The electric current density for |$sDw) is calculated as 



j(r) = -i^— (isDV^|^i(r)VV'a(r)|4>SDw) 

Hp 



m 

a 



. Vx, (18) 
2m 

where m is an effective mass, p is the particle density, and the current from the ordinary SDW 
state without merons is calculated to be zero. 

Let us examine a current generation by a collection of merons, 

-1 y-Vc 

c X Xp. 

which gives rise to spin density modulation. 



e(r) = ^tan-i^, (19) 



Sxir) = ($sDw|^ (v5{(r)V'i(r) +V'I(r)V^T(r)) I^SDw) 



'2 

= |s|cos(QT + ^(r)), (20) 
Sy{v) = ($sDwl^ (^^}(r)^^(r) - '0|(r)^^(r)) |$sDw) 

= |s|sin(QT + ^(r)), (21) 
s^(r) = (l-sDwl^ (^}(r)V;|(r) - V;{(r)^^(r)) |$sDw) 

= 0. (22) 
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Fig. 2. Top: a spin texture generated by three merons and an SDW order with Q — (tt, tt) in the 
two-dimensional space. The left solid line indicates the left edge of the system. Bottom: a current 
induced by merons, where x is chosen to be x = 

An example of a two dimensional spin texture generated by an SDW order with three 
merons is depicted in Fig. 2. A corresponding current is also shown in Fig. 2 by taking x = ^• 

If the balance of upward and downward currents is destroyed, a collection of meron currents 
may give rise to a macroscopic current. To investigate such a possibility, we first examine the 
meaning of the nontrivial phase x- 

An appearance of the nontrivial phase actually means an appearance of a magnetic field. 
The phase factor e^a^^'") changes the derivative to a covariant derivative, 

V/(r)^ (^V + iivx)/(r), (23) 

which may be compared with the covariant derivative in the presence of a magnetic field 
B = V X A, 

(V-,^a)/M. (24) 

If we define a fictitious vector potential Age by 

Afie = -^Vx, (25) 
zq 

the correspondence between the r.h.s of Eq. (23) and Eq. (24) suggests that. Age plays a role 
of the vector potential for a magnetic field. 

It is known that the effect of a magnetic field on wave functions amounts to assign a non- 
integrable phase factor to every loop along which the charged particle travels in semiclassical 
sense. The information about this is encoded in A as a connection. The vector potential 
Afic does the same assignment for electrons in the SDW band. 
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In other words, for the SDW electronic state with merons that give rise to the phase ^, 
the wave function is the one augmented by a nonintegrable phase factor that is exactly the 
same one produced by the magnetic field given by V x Age- Thus, as to the effect on electrons 
in the SDW band, Age plays the same role as a true vector potential, and V x Age as a true 
magnetic field. 

Therefore, we have 

he 

B = V X Afic = -—V X Vx, (26) 
Iq 

which is not zero due to the multi-valuedness of %. It is also worth mentioning that if Age is 
given by — yVx instead of — |^Vx, the phase factor for every loop is +1 (the phase is 27rn 
where n is an integer); thus, a nonintegrable phase factor does not arise. 
Using Afic, Eq. (18) is written as 

j(r) = -^Afie. (27) 
mc 

This is the London equation^^ with the real vector potential replaced by the fictitious one. 
Prom Eqs. (26) and (27) , and one of Maxwell's equations 

V X B = ^j, (28) 

we obtain an equation derived by London:^^ 

V'B(r) = ^B, (29) 

where Afic(r), B(r) and j(r) denote averages of Afic(r), B(r) and j(r) over a small region 
around r (we will denote the similar average value for a quantity Q as Q), and A/, is the 
London penetration depth^^ given by 



As is well-known Eq. (29) describes an expulsion of a magnetic field. 

An equation that relates the average magnetic field to the average meron current is given 
from Eq. (26) as 



B = X Vx, (31) 



where $o is a fluxoid given by 



*. = |. (32) 



Here, electron charge — e = q\s used. 

Equation (31) has a simple meaning: if x Vx is integrated over a surface S with a 
boundary C, we have 

i/^VxVx.dS = i£vx-<ir = i^£<iX, (33) 
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Fig. 3. A current distribution generated by sheets of meron current centers. The system is uniform 
in the z direction. Each sheet of meron current centers is represented by a dotted Hne. It causes a 
jump in the current. The current is a diamagnetic one, which decays to zero as the distance from 
the surface increases. 



where the last expression is equal to the difference of the number of counterclockwise meron 
current centers (each center contributes +1) and that of clockwise ones (each center contributes 
— 1) within the loop C. We call the quantity in Eq. (33), " the sum of meron current winding 
numbers (SMCW)". Therefore, Eq. (31) is actually ^>o times an average of the SMCW in 
units area normal to the magnetic field. 

Equations (28), (29), and (31) comprise a system of equations for a stable macroscopic 
current. Note that when the three equations are consistently satisfied, B is considered to be 
a true magnetic field produced by j through Eq. (28). 

Let us obtain a solution for a semi infinite system that extends in the x > region 
(Fig. 3). A magnetic field points in the z direction. The average magnetic field is calculated 
from Eq. (29) as 

B,(x) = S,(0)e"^, (34) 

which decays as x increases. 

Prom Eq. (31), an average of the SMCW is obtained as 

Pm(x) = ^e-^. (35) 

According to Eq. (28), the diamagnetic current flows in the y direction, and its relation 
to pm is obtained as 

The diamagnetic current may be considered to be generated by sheets of meron current 
centers: from Eqs. (25) and (27), we have 



Vxj = -^3^Vx Vx, (37) 
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which is analogous to an equation in hydrodynamics given by 

V X u = a;, (38) 

where u and uj are velocity and vorticity in a fluid, respectively. We see correspondence 
between j and u, and V x Vx and uj. It is known in hydrodynamics that a vortex sheet 
causes a jump in the fluid velocity/^ which in the present case means that a sheet of meron 
current centers causes a jump in the current density. 

The current should be zero deep inside the sample; a small current that would flow if the 
electric resistivity of the sample is zero does not arise in real materials. Starting from there, a 
successive positive current jumps through sheets of meron current centers toward the surface 
leads to a current density distribution that decays as x increases (Fig. 3). Actually, the current 
density jump is an extreme ideal case; in reality it varies gradually and continuously. 

So far we have been considering the case with a collection of merons only; however, the 
formation of meron-antimeron pairs is energetically more favorable since the background SDW 
order is not destroyed in the region away from them.^'^ Each meron current around a meron 
or an antimeron has two choices in its directions. When a macroscopic current is present, 
one of them is chosen. Merons or antimerons are not easily destroyed at sufficiently low 
temperature because they are topological objects, thus the destruction of them requires a 
cooperative movement of a large number of electrons. Major dissipative processes for the 
diamagnetic current will be single particle excitations and meron flow. When these processes 
are negligible, the diamagnetic current should persist. 

In summary, we have investigated a system of an SDW order with merons. It is shown 
that a stable circular current appears around a meron. A collection of such currents will gives 
rise to a diamagnetic current. If the meron density is sufficiently large, a perfect diamagnetism 
is realized. The stability of the diamagnetic current is attributed to the topological nature 
of merons, and if dissipative processes such as single particle excitations and meron flow are 
negligible, the current will be a persistent one. 
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